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Abstract 


Minimization  of  free  energy  is  used  to  calculate  the  equilibrium  shapes  of  sessile  liquid 
drops  under  different  wetting  conditions  and  under  the  influence  of  a gravitational 
body  force.  The  total  system  free  energy  is  assumed  to  obey  an  equation  recently 
derived  by  Searcy  [1],  and  the  shapes  of  sessile  drops  of  different  sizes  are  calculated 
and  compared  with  experimental  measurements.  Stationary  states  in  the  free  energy  are 
found  using  standard  variational  principles,  which  lead  not  only  to  an  Euler-Lagrange 
differential  equation  for  the  drop  shape,  but  also  to  the  boundary  conditions  at  the 
three-phase  junction  where  the  liquid,  solid,  and  surrounding  vapor  meet.  The  analysis 
shows  that  the  Searcy  equation  leads  to  the  classical  Young-Dupre  equation  for  the 
thermodynamic  contact  angle.  Drop  shapes  calculated  by  this  approach  are  in  excellent 
agreement  with  experimental  observations  of  drops  of  pure  water  resting  on  a clean 
polytetrafluoroethylene  (PTFE)  substrate. 
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1 Introduction 


Wetting  and  adhesion  phenomena  are  the  basis  of  a diverse  range  of  natural  phenomena  and 
engineering  processes.  The  tendency  of  a liquid  to  wet  a solid  is  determined  by  the  free  energy 
reduction  that  accompanies  the  covering  of  the  solid  by  the  liquid,  relative  to  the  combined  free 
energies  of  the  same  mass  of  liquid  and  solid  in  isolation. 

A large  number  of  theoretical  and  experimental  investigations  of  wetting  have  been  published 
in  the  past  200  years.  Many  of  these  have  focused  on  the  equilibrium  shape  of  a liquid  drop  or 
a solid  particle  resting  atop  a planar  solid  surface.  A schematic  of  a sessile  liquid  drop  partially 
wetting  a flat  solid  is  shown  in  Fig.  1(a).  If  the  drop  is  sufficiently  small,  or  if  the  influence  of 
gravity  is  negligible,  then  the  drop  is  observed  to  have  a shape  that  approximates  a spherical  cap, 
although  the  particular  shape  is  not  of  concern  at  this  point.  What  is  important  is  that  the  angle, 
Ow,  at  which  the  liquid- vapor  surface  meets  the  solid-vapor  surface  appears  to  be  a thermodynamic 
property  depending  only  on  the  three  different  interfaces  that  are  present  (liquid- vapor,  solid- vapor, 
and  liquid-solid).  This  angle  often  is  called  the  thermodynamic  contact  angle  for  the  system. 


(a) 


Figure  1:  (a)  Schematic  drawing  of  a sessile  liquid  drop  partially  wetting  a flat  solid  substrate  with 
thermodynamic  contact  angle  0^.  cj  indicates  the  gravitational  acceleration  vector,  (b)  Axisym- 
metric  liquid  drop,  showing  the  spherical  polar  coordinate  system  that  is  used  in  the  derivation. 

The  accepted  relation  between  the  thermodynamic  contact  angle  and  the  thermodynamic  prop- 
erties of  the  surfaces  is  the  Young-Dupre  equation: 


^ 7 SV  lls  / "i  \ 

cos  0^  = (1) 

llv 

where  gsv,  q/s,  and  giv  are  the  surface  free  energy  densities  of  the  solid- vapor  surface,  the  liquid-solid 
surface,  and  the  liquid-vapor  surface,  respectively. 

At  constant  temperature  and  volumes  of  the  three  phases  in  Fig.  1(a),  the  thermodynamic 
potential,  G,  which  governs  equilibrium  and  stability  of  the  system,  is  [3] 


G — Gjj  + 7 i-A-i 

i 
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where  Gj,  is  the  contribution  to  G of  the  bulk  phases,  Ai,  is  the  area  of  the  i-th  type  of  surface,  and 


li  = 


Aj  ^A, 


is  the  surface  energy  density  of  the  z-th  type  of  surface.  The  summation  in  Eq.  (2)  represents  the 
contribution  to  the  thermodynamic  potential  of  the  excess  free  energy  stored  in  the  surfaces. 

At  thermodynamic  equilibrium  any  inhnitesima]  change  in  the  configuration  of  the  system 
cannot  decrease  G.  For  changes  at  constant  temperature  and  constant  volume  of  the  liquid,  solid, 
and  vapor,  G 5 is  constant  and  the  condition  of  equilibrium  is 


SG  — S (7 ivAiv  + 'jsyAsy  + 7 isAis)  > 0 (3) 

which  is  a particular  form  of  the  more  general  Gibbs  criterion  [3,4]. 

Recently,  Searcy  has  proposed  an  alternative  to  Eq.  (3)  for  the  condition  of  equilibrium, 


5G{  — 5 [7 ivAiv  + (7 iv  + A 1)  Aig]  > 0 (4) 

where  now  the  potential  G'/  is  the  free  energy  of  the  liquid,  and  A 1 is  the  partial  free  energy  (energy 
per  mole  of  liquid)  of  bonding  of  the  liquid  to  the  solid,  per  unit  area  of  the  liquid-solid  surface 
that  is  formed.  The  parameter  A 1 can  be  interpreted  in  the  following  way.  Imagine  that  a liquid 
drop  is  somehow  deformed  into  exactly  the  shape  it  will  have  when  it  wets  the  substrate,  but  that 
it  is  not  yet  brought  into  contact  with  that  substrate.  The  free  energy  of  the  drop  in  this  deformed 
state,  relative  to  the  bulk  liquid,  is  well-defined  and  depends  only  on  its  surface  area  and  7 iv.  Now 
the  deformed  drop  is  brought  into  contact  with  the  solid  substrate.  The  shape  of  the  drop  does  not 
change,  so  the  only  change  in  free  energy  is  due  to  the  chemical  bonding  of  liquid  and  solid  across 
the  newly  formed  surface.  The  parameter  A 1 is  equal  to  the  difference  between  the  free  energy  of 
the  drop  in  the  latter  case  and  the  free  energy  of  the  identically  deformed  drop  without  wetting, 
per  unit  area  of  contact  between  liquid  and  solid.  The  free  energy  change  is  assigned  solely  to  the 
liquid,  while  the  free  energy  of  the  substrate  is  assumed  to  remain  constant  . 

Both  Eq.  (3)  and  Eq.  (4)  ignore  the  influence  of  gravity  on  the  thermodynamics  of  the  system, 
and  it  is  well  known  that  the  equilibrium  shape  of  a macroscopic  sessile  drop  is  sensibly  affected 
by  gravity.  Most  notably,  the  equilibrium  shape,  which  is  size-invariant  in  the  absence  of  gravity, 
is  dependent  on  the  volume  of  the  drop  when  gravity  is  present. 

The  purpose  of  this  paper  is  to  provide  a thermodynamic  basis  for  calculating  the  equilibrium 
shape  of  an  axisymmetric  sessile  liquid  drop  under  the  influence  of  gravity.  That  is,  the  approach 
used  here  relies  011  a formal  minimization  of  the  appropriate  thermodynamic  potential,  using  varia- 
tional calculus.  The  analysis  leads  to  an  Euler-Lagrange  differential  equation  describing  the  shape 
of  the  drop.  In  addition,  the  approach  also  naturally  derives  the  boundary  condition  that  must 
prevail  at  the  three-phase  junction.  Prior  investigations  similar  to  this  have  usually  assumed  the 
boundary  condition  according  to  Eq.  (1).  Therefore,  the  approach  used  here  seems  more  satisfying 
from  a thermodynamic  perspective.  Furthermore,  this  approach  provides  a means  to  directly  test 
the  form  of  the  thermodynamic  potential  proposed  by  Searcy  in  Eq.  (4).  It  will  be  shown  that  the 
implications  of  Searcy’s  equation,  at  least  for  the  shape  and  properties  of  macroscopic  sessile  drops, 
is  indistinguishable  from  those  of  the  accepted  thermodynamic  potential  originally  formulated  by 
Gibbs,  a special  case  of  which  is  given  in  Eq.  (2). 
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2 Theory 

2.1  Derivation  of  Governing  Equations 

The  system  being  considered  is  depicted  in  Fig.  1(a).  It  consists  of  a liquid  drop  having  constant 
volume  V immersed  in  a vapor  phase  and  resting  on  a planar  solid  substrate.  The  drop  is  assumed 
to  have  an  axis  of  rotational  symmetry  represented  by  the  vertical  dashed  line  in  Fig.  1(b).  A 
gravitational  body  force  g is  assumed  to  act  in  the  —xs  direction. 

For  purposes  of  being  explicit,  the  following  assumptions  are  made  about  the  system  before 
proceeding  with  the  analysis. 

Assumption  1 The  solid  is  rigid  and  has  a planar  interface  with  the  vapor  and  liquid. 


Assumption  2 Each  phase  has  fixed  composition  and  can  undergo  no  chemical  reactions  or  segre- 
gation. 


Assumption  3 The  surface  energy  density  of  each  surface  is  constant.  In  particular,  there  is  no 
crystalline  anisotropy  in  the  surface  energy  density  of  any  of  the  interfaces. 


Assumption  4 Each  phase  has  constant  mass  and  volume. 

This  assumption  is  reasonable  for  the  kinds  of  wetting  phenomena  that  are  being  examined  here. 
It  should  be  mentioned  that  this  assumption  prohibits  the  formation  or  annihilation  of  vacancies 
in  the  solid  or  liquid  that  can  change  the  volume  at  constant  mass.  However,  the  contribution  to 
the  configurational  component  of  the  thermodynamic  potential  that  these  process  would  allow  is 
usually  small  and  insensitive  to  changes  in  shape  at  constant  mass. 


Assumption  5 The  vapor  is  a thermal  reservoir  that  fixes  the  temperature  of  the  liquid  and  solid 
phases. 

Assumption  6 The  liquid  drop  is  rotationally  symmetric  about  an  axis  which  intersects  the  centroid 
of  the  drop. 

Assumption  7 The  gravitational  field  is  constant  over  the  volume  of  the  drop. 

The  thermodynamic  potential  is  composed  of  a contribution  G5  from  the  bulk,  a net  contribution 
Gi  from  the  interfaces,  and  a contribution  Gg  due  to  gravitational  potential  energy: 


G — Gb  + Gi  + Gg 

For  the  system  in  Fig.  1,  Gi  can  be  written  according  to  Gibbs  [3], 


(5) 
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G%  'llv-A-iv  “I-  ^fsv-^-sv  T lls-^ls 


(6) 


As  stated  in  the  Introduction,  Searcy  [1]  proposed  the  following  modified  equation  for  the  system 
in  Fig.  1 


Gi  — 7 ivAiv  + (7 iv  + ^1)  Ais  (7) 

For  the  remainder  of  this  section,  we  will  use  Eq.  (7)  to  represent  Gi , and  will  compare  the  results 
to  those  that  can  be  obtained  when  Eq.  (2)  is  used  instead. 

The  contribution  to  G of  gravitational  potential  energy  is  given  by  an  integral  over  the  volume 
of  the  drop: 


Gg  = - J pg  •f{xi,X2,xs)dV  (8) 

where  r is  the  position  vector  relative  to  the  origin  O,  g is  the  acceleration  due  to  gravity,  and  p is 
the  density  of  the  liquid.  Assumptions  4 and  7 allow  us  to  simplify  Eq.  (8)  to 


Gg  = pg  x3  ■ r{x1,x2,x3)dV  (9) 

Jv 

where  now  g = \g\  is  the  magnitude  of  the  acceleration  due  to  gravity.  The  integral  in  Eq.  (9)  will 
be  recognized  as  the  expression  for  the  elevation  of  the  centroid  of  the  drop  above  the  plane  of  the 
substrate. 

Because  of  Assumption  1,  we  may  conveniently  use  a spherical  polar  coordinate  system,  as 
shown  in  Fig.  1(b).  The  rotational  symmetry  of  the  drop  means  that  its  shape  is  independent  of 
the  azimuthal  angle  (J).  Thus  the  liquid-vapor  interface  may  parametrized  by  the  polar  angle  9.  so 
that  the  ordered  pair  of  parametric  functions  {X\  (#),  X3(9))  specifies  the  location  of  liquid- vapor 
interface  for  any  value  of  9.  X\  and  X3  both  are  one-to-one  functions  of  9 unless  some  portion  of  the 
liquid- vapor  interface  lies  parallel  to  the  position  vector  R{9)  in  Fig.  1(b)  for  9 G (0,7r/2).  Later, 
in  Section  3,  a different  parameterization  will  be  helpful  for  numerical  solution  of  the  problem. 
However,  for  now  the  parameterization  in  9 is  convenient  for  the  variational  derivation  because  it 
provides  fixed  minimum  and  maximum  values  of  the  parameter  (0  and  7t/2,  respectively). 

With  the  coordinates  system  defined,  for  now,  as  in  Fig.  1(b),  we  can  rewrite  Eq.  (7)  as 

Gi  = 7r  (7i„  + A,)  Xl  (|)  + 2nllv  G"  X^e) /wT)  + AT)  d 9 (10) 

where  X\  and  X3  represent  the  first  derivatives  of  X[(0)  and  X3(9),  respectively.  Also,  X\  and  X3 
both  can  be  written  in  terms  of  R in  Fig.  1(b): 


X! 

= RsinO 

X\  = R sin  9 + R cos  6 

(ii) 

X3 

= R cos  6 

X3  = R cos  9 — R sin  9 

(12) 
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Variable  Dimensionless 

Variable  Dimensionless 

A,  a = ^~ 

1 7 lv 

II 

& 

G Q=  G, 2 

7T7  IvL2 

V V=  ^3 

I 7T/V 

Table  1:  Dimensionless  groupings  of  variables. 


Therefore,  Eq.  (10)  becomes 


Gt  = w (7|„  + A,)  R2  (|)  + 2rr7i„  R(6) \j  R2(0)  + R2(0)  sin 9 d W (13) 


The  gravitational  term  Gg  in  Eq.  (9)  can  be  explicitly  written  in  integral  form  as 


rir/2  pR(0) 

Gg  — 2irpg  / sin#  d#  / (rcos#)r2dr 

Jo  Jo 

rn/2 

— ~P9  / F4(#)  cos# sin#  d#  (14) 

2 Jo 


The  thermodynamic  potential  in  Eq.  (5)  can  now  be  written  in  integral  form  as 


G = g(r,R}  =ir('yiv  + \i)R2{tt/2) + ir  j [2 ~nvRpR2  + R2  + ^ gR 4 cos sin0  dd  (15) 


We  seek  to  minimize  F subject  to  the  constraint  that  the  volume  V of  the  liquid  is  fixed: 


L3  = V 


2tt 

T 


•*7t/2 


F3  sin  # d#  = constant 


(16) 


where  we  dehne  L — E1/3  to  be  a characteristic  dimension  of  the  drop. 

The  analysis  is  expedited  by  defining  dimensionless  variables  as  shown  in  Table  2.1.  With  these 
dimensionless  groupings,  Eqs.  (15)  and  (16)  become 


(71.71)  = {l  + a)K~  (- ) + 


7T 


’7r/2 


2fVf2  + F2  + ^F4  cos 


sin#  d# 


(17) 


2 W2 

V = — / F3sin#d#  (18) 

3 Jo 

The  variational  approach  used  here  is  adapted  from  Segel  [5].  Assume  that  F(#)  is  the  extrem- 
izing  function,  and  examine  small  perturbations  of  this  function  having  the  form 
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K(e)  + eMO)  + e&m 


(19) 


51(0) 

51(0) 


7t(0)  + €1C1(0)  + e2C2(0) 


(20) 


where  (j  and  C,2  are  arbitrary  C\  smooth  functions  of  9 and  ei,  e2  <C  1. 

We  can  substitute  the  trial  solution,  Eqs.  (19)  and  (20)  into  Eqs.  (17)  and  (18)  to  obtain 


g{e1,e2)  = (l+a)tf 


(!)+. 

r 7r/2 

lo 

291\Ak2  + 912  + 

— 914  cos  9 
2 

sin  9 d 9 

(21) 

2 

r/2  q 

(<T,  ^2) 

~*J 

913  sin  9 d 9 

0 

(22) 

Using  the  method  of  Lagrange  undetermined  multipliers,  the  minimzation  of  Eq.  (21)  subject  to 
the  constraint  in  Eq.  (22)  is  equivalent  to  the  minimzation  of  a new  function: 


r (61,62)  = Q (61,62)  - V V(ei,e2) 


(23) 


where  i)  is  an  undetermined  constant.  The  function  T is  a function  only  of  the  fluctuations  e\  and 
e 2 . A stationary  value  of  T must  satisfy 


ar 

dei 


£2 


ar 

de2 


ei 


For  convenience  in  subsequent  equation,  we  adopt  the  shorthand  notation 


Ti  = (llWlH2  + 5t2  + 1<K4  cos  0^  sin  0 


V2  = -943sim 
3 


Then,  evaluating  the  left  derivative  in  Eq.  (24), 


(24) 


(25) 

(26) 


ar 

dci 


->  ei,e2=0 


-2a+„,«,|)=iil 


+ 


ei,£2=0 

fn/2/'drl  dr2 

0 


de-i 


V 


dci 


&9 


0 z = l,  2 (27) 


J ei,e2=0 


The  derivatives  inside  the  integral  can  be  expanded  by  the  chain  rule, 
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j = 1 » 2;  z — 1,2 


(28) 


ari_ar1cm  ar.dSR  dr?  d# 

dej  <991  dfj  <991  dej  90  de3 

Because  0 is  the  independent  variable,  dO/dej  = 0,  so  the  third  term  in  this  equation  is  zero. 
At  the  endpoints  0 = 0 and  0 — tt/2,  the  91  varies  with  e\  and  e2  according  to 


d9t(0) 

dei 

d91(0) 

de2 

dW(f) 

dei 

dW(f) 

de2 

In  addition,  Eq.  (19)  implies  that 


Ci  (o) 

(29) 

(2(0) 

(30) 

to|  ^ 

(31) 

& (?) 

(32) 

91  ( - 


9 

d£i 

<991 


= ni- 


ei,e2=0 


ei,e2=0 


dT \ 

dn  ’ 


and  finally 


di\  =d£i 

d% 1 d ,e2=0  dJl 


2 — 1,2 


(33) 


Substituting  Eqs.  (29)-(33)  into  Eqs.  (23),  (24)  and  (28)  gives 


dT 

dej 


= 0 = 2(l  + a)ft  - Ki  - + 


7T 


7T 


^ r/^ 


ard 

dn 


did 

' dn 


Ci  d 0 


-tt/2 


+ 


di\ 

dn 


or  2 

' dh 


CidO  (34) 


The  second  integral  can  be  rewritten  using  integration  by  parts: 


-tt/2 


dli 

dn 


dr  2 
dn 


/ 7T  \ 

Ci  do  = 

V 2 / 

dT 1 

dn 

C<(0) 


6>r.; 


0- 


dn 


- 6=tt/2 


<9r2l 

j-TT/ 

2 d 

r^Ti 

dry 

.art 

/;  <97?.. 

Q} 

II 

O 

1 

O " 

dd 

,<97t 

/ / 

d7?J 

Ci  d 0 (35) 


Also,  from  Eq.  (26),  it  is  clear  that  dT2/On  = 0.  Therefore,  substituting  and  grouping  like  terms 
in  Ecp  (34)  gives 


/ 7T  \ 

, . „ /7T\  dTl 

2 1+a  R - +4 

-C:(0) 

[^1 

- 

V 2 / 

L V2/  on \ 

9=ir/2 

[on  \ 

0T1  d /dTih 
dU  d 0\dil) 


dY, 


li- 


on 


Ci  m 


2 = 1,2  (36) 


Because  Ci  and  C2  are  arbitrary  smooth  functions,  the  only  way  to  guarantee  that  Eq.  (36)  vanishes 
is  if  each  term  separately  is  zero: 


. / it  \ 5r  1 

2 (1  + ck)  ft  ( — ) + — — 

v 2 / on 


6=ir/2 


(37) 


or 1 
on 


6= 0 


= 0 


(38) 


»7r/2 


r^r! 

d 

40  \ 

or  2 1 

on 

d e 

Koh) 

''a  ft 

Cl  d6»  = 0 


E 2 


(39) 


Ecp  (39)  was  derived  under  the  assumption  that  Ci  and  C2 , while  arbitrary,  are  constant.  Therefore, 
in  principle,  77  could  depend  on  them.  However,  it  may  be  shown  [5]  that  //  must  be  independent  of 
Ci  or  C2  or  both.  Therefore,  the  term  in  square  brackets  of  the  integrand  in  Ecp  (39)  must  be  zero: 


0Tl  d /<9Ti\  0Y2 

~&R~dd  V 97t  ) ~ 


(40) 


To  make  the  notation  more  compact,  let 


/ (ft,  ft)  = r,  — r 2 

Then  Eqs.  (37),  (38),  and  (40)  take  the  form 


(41) 
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A (K 
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= 0 
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(42) 

(43) 


(44) 


2.2  Evaluation  of  Governing  Equations 

Before  solving  Eq.  (44),  we  can  gain  insight  about  the  form  the  solution  must  have  just  by  further 
evaluating  the  two  boundary  conditions,  Eqs.  (42)  and  (43),  and  the  Euler-Lagrange  equation  (44), 
in  terms  of  7 Z.  In  particular  we  will  see  that  the  boundary  conditions  specify  the  thermodynamic 
contact  angle  at  the  three-phase  junction  ( 0 — tt/2) , and  that  the  Euler-Lagrange  equation  provides 
a physical  meaning  to  the  Lagrange  multiplier,  77. 


2.2.1  Evaluation  of  Boundary  Condition  at  9 = 0 

Substituting  for  / in  Eq.  (42)  gives 


df 

21Z1Z 

dn 

d=0 

yj  1Z2  + 772 

2ft(0)7^(0) 

sjii2(0)  + n2(<a) 


(45) 


From  Fig.  1(b)  we  require  that  77.(0)  > 0 for  finite  drop  volumes.  Therefore,  Eq.  (45)  implies  that 


7t(0)  = 0 Boundary  Condition  1 (46) 

which  means  that  the  drop  profile  is  flat  at  its  apex.  This  is  a natural  consequence  of  the  rotational 
symmetry  that  has  been  imposed  in  the  statement  of  the  problem.  Under  rotational  symmetry, 
any  nonzero  value  of  1Z  at  6 = 0 would  correspond  to  a singularity  in  mean  curvature  which  is 
unphysical,  as  we  will  see  later. 

2.2.2  Evaluation  of  Boundary  Condition  at  6 = tt/2 
Substituting  for  / in  Eq.  (43)  gives 


21 

dii 


B=tt/2 


2TC(f)7t  (f) 

/v2  (I)  + V2  (f) 


-2(1  + a)Tl{^) 


which  can  be  simplified  to 


+ 7^2 


Squaring  both  sides  and  rearranging  leads  to  a further  simplification 


• / 7T  \ (1  + (| 

n(  -)  = ± , 

'2/  v/-(2a  + a2) 

But  in  Eq.  (48),  the  radical  on  the  right  must  be  non- negative,  so  we  expect  that 

1l  < 0 a > -1 

Therefore,  we  choose  the  negative  root  in  Eq.  (49). 


(47) 


(48) 


(49) 
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Figure  2:  Enlarged  view  of  region  near  the  three-phase  junction  between  solid,  liquid,  and  vapor. 


h 


(1  + a)  TZ  (f ) 
\f  (2 ot.  -f-  cv" ) 


(50) 


This  equation  places  limits  on  the  possible  values  of  a,  because  TZ(7t/2)  is  required  to  be  a real 
value.  Therefore,  a condition  on  a is 


— 2 < <a  < 0 


or,  equivalently 


-2ltv  < \i  < 0 


(51) 


(52) 


Furthermore,  1Z(t:/2)  can  be  related  to  the  thermodynamic  contact  angle,  0w;,  in  Fig.  1(a).  An 
enlarged  view  of  the  region  near  the  three-phase  junction  is  shown  in  Fig.  2.  referring  to  the  small 
triangle  outside  the  drop  in  this  figure,  basic  trigonometry  requires  that 


U ( 9 ) cos  (A 9)  -H{6  + A6)  = Vil2  + 7Z2  cos  (tt  - Qw)  A 6 (53) 

Dividing  both  sides  of  Eq.  (53)  by  A0,  taking  the  limit  as  A 9 approaches  zero,  and  evaluating 
at  0 = 7r/2, 


1 im  \/  R 2 + 'R  - cos  ( 7r 
Ad^o 


07 


- 6=ix/2 


lim 


1Z{9)  -K{0  + AO) 


a e 


- 0= 7r/2 


-n 


Therefore,  using  the  fact  that  cos(7r  — x)  = — cos(j>),  we  have 


MS) 

(f ) + K2  (f ) 


cos  Qw 


(54) 


(55) 
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Comparing  to  Eq.  (48)  shows  that 


cos  ©u;  = — (1  + a)  Boundary  Condition  2 (56) 

Substituting  for  a from  Table  2.1  gives: 


liv  cos  — llv  T A i (57) 

Searcy  gives  this  equation  in  the  same  paper  in  which  he  introduced  the  partial  free  energy  of 
bonding,  A i [1],  and  it  now  has  been  shown  how  it  arises  naturally  as  a boundary  condition  on  the 
solution  for  the  free  energy-minimizing  drop  shape. 


Remark  1 If  we  make  the  assignment 


A i — ('7  is  Ksv)  llv  (58) 

then  Eq.  (57)  is  identical  to  the  Young-Dupre  result  in  Eq.  (1). 


Remark  2 Eq.  (57)  is  independent  of  the  influence  of  gravity,  even  though  gravitational  contribu- 
tions to  the  thermodynamic  potential  have  been  included  in  the  problem. 

2.2.3  Evaluation  of  the  Euler- Lagrange  Equation 

Substituting  for  / in  Eq.  (44)  gives,  after  a considerable  amount  of  algebraic  manipulation, 


2 n2  + ft2  -iin  n-n  cot  e 

{fr  + wf2  + nViv  + iv 


+ (3 1Z  cos  6 — p 


(59) 


Recall  that  r/,  the  Lagrange  multiplier  imposed  by  the  volume  constraint  on  the  drop,  must  be 
constant.  Therefore  the  grouping  of  terms  on  the  left  side  of  Eq.  (59)  must  have  fundamental 
significance.  In  fact,  hZ{9)  cos  0 is  the  dimensionless  elevation  (h/L)  above  the  substrate  surface  at 
a point  on  the  drop  surface  specified  by  6.  And  because  in  Table  2.1  (3  is  the  dimensionless  form 
of  pg,  the  last  term  on  the  left  side  of  Eq.  (59)  is  the  dimensionless  contribution  of  the  hydrostatic 
pressure  at  any  height  within  the  drop  due  to  gravity. 

Principles  of  differential  geometry  [6]  can  be  used  to  show  that  the  term  within  square  brackets 
on  the  left  side  of  Eq.  (59)  is  the  local  (dimensionless)  mean  curvature  of  the  liquid-vapor  sur- 
face1, K,  = LK,  at  any  point2  specified  by  6 (see  Appendix  A).  Therefore,  this  term  represents  the 
contribution  to  the  pressure  within  the  drop  due  to  the  presence  of  a curved  surface  as  originally 
derived  by  Laplace  [7]. 

In  summary,  then, the  form  of  the  Euler-Lagrange  equation  is  equivalent  to 

1 We  use  the  convention,  common  in  the  theory  of  capillary  phenomena,  of  denoting  by  the  term  “mean  curvature” 
a property  of  surfaces  that  is  actually  two  times  the  mean  curvature  as  it  is  defined  in  differential  geometry. 

2For  0 — > 0,  cot#  in  Eq.  (59)  is  singular.  Appendix  A deals  with  this  in  greater  detail. 
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(60) 


At  any  point  on  the  liquid-vapor  surface,  then,  the  mean  curvature  must  decrease  as  the  height 
increases  in  such  a way  as  to  keep  the  pressure  within  the  drop  constant  everywhere.  Thus,  77  is  seen 
to  be  the  thermodynamic  force  conjugate  to  the  volume  of  the  liquid,  which  is  constant  everywhere 
inside  the  drop  at  mechanical  equilibrium. 

Eq.  (60)  was  derived  by  minimizing  Searcy’s  expression  for  the  free  energy  (Eq.  (7).  However, 
because  Eq.  (7)  and  Eq.  (2)  have  exactly  the  same  functional  dependence  on  the  areas  of  the  various 
interfaces,  Eq.  (60)  does  not  depend  on  which  of  these  two  expressions  for  the  thermodynamic 
potential  is  used. 

3 Numerical  Results 

The  remaining  task  is  to  solve  the  Euler-Lagrange  equation  (59)  subject  to  the  boundary  conditions 
given  in  Eqs.  (45)  and  (50).  Eq.  (59)  is  a second-order  nonlinear  ordinary  differential  equation 
(ODE),  and  there  are  no  analytical  techniques  available  to  solve  it.  Therefore,  we  appeal  to  a 
numerical  solution.  To  numerically  solve  an  nth-order  ODE,  it  is  first  convenient  to  re-express  it 
as  n coupled  first-order  ODEs,  and  there  are  standard  techniques  to  accomplish  this.  Furthermore, 
Eq.  (59)  has  a singularity  at  0 = 0.  The  singularity  can  be  handled  by  reparameterizing  the 
problem  in  terms  of  arc  length  of  the  liquid-vapor  surface  and  by  transforming  from  spherical  polar 
coordinates  to  cylindrical  coordinates.  When  this  is  done,  the  singularity  is  still  present  but  can  be 
removed  by  expanding  about  the  singular  point.  The  details  of  the  reparameterization,  coordinate 
transformation,  and  removal  of  the  singularity  are  given  in  Appendix  B. 

The  numerical  approach  used  here  involves  two  algorithms  from  the  NIST  core  mathematical 
library,3  (“CMLIB”):  SDRIV3  and  SNSQE.  These  use  a shooting  method  [8]  from  the  apex  of  the 
drop  down  to  the  contact  line  with  the  substrate.  For  a given  drop  volume,  the  solution  is  iterated 
by  adjusting  the  value  of  the  Lagrange  multiplier  until  the  desired  drop  volume  is  obtained. 

In  the  results  that  follow,  we  relate  the  equilibrium  shapes  of  sessile  drops  to  physical  quantities 
instead  of  the  dimensionless  parameters  used  in  the  previous  section.  To  that  end,  we  assume  values 
°f  llv  — 72.2  mJ/m2  and  p = 1000  kg/m3,  which  are  consistent  with  the  value  reported  for  pure 
water  in  contact  with  air.  We  also  assume  g = 9.8  111/s2. 

Figure  3 shows  calculated  cross-section  profiles  of  a sessile  liquid  drop  with  volume  of  1 mL.  The 
profile  is  shown  for  three  different  values  of  the  thermodynamic  contact  angle  . The  main  effect 
of  increasing  Qw  is  to  pull  the  drop  into  a more  equiaxed  shape  having  greater  mean  curvature. 

The  influence  of  gravity  on  the  equilibrium  shape  becomes  increasingly  pronounced  as  the 
drop  volume  increases.  Fig.  4 shows  the  calculated  equilibrium  profiles  for  different  drop  volumes 
when  Qw  = 90°.  Larger  drops  adopt  more  oblate  shapes  relative  to  smaller  drops  because  the 
gravitational  contribution  to  the  free  energy  is  reduced  by  having  the  liquid  distributed  at  lower 
elevations  above  the  substrate. 

It  is  interesting  to  compare  the  calculated  equilibrium  profiles  of  drops  to  experimental  obser- 
vations of  sessile  drops.  Recently,  Barberis  and  Beruto  performed  measurements  of  water  drops 
resting  on  a flat  surface  of  polytetrafluoroethylene  (PTFE)  [2].  Both  the  water  and  PTFE  were 
carefully  treated  to  avoid  surface  contamination.  Measurements  made  of  the  water-air  surface 

'The  NIST  core  mathematical  library  can  be  found  at  the  NIST  Guide  to  Available  Mathematical  Software, 
http : / /gams . nist . gov/. 
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Figure  3:  Cross-sectional  profiles  of  equilibrium  shapes  of  liquid  sessile  drops  with  volume  of  1.0  mL. 
Profiles  are  shown  for  Qw  = 45°,  90°,  and  135°. 


Figure  4:  Cross-sectional  profiles  of  equilibrium  shapes  of  liquid  sessile  drops  with  0^  = 90°. 
Profiles  are  shown  drop  volumes  of  0.1  mL,  0.5  mL.  1.0  mL,  and  2.0  mL. 
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Figure  5:  Photograph  of  a 2 niL  drop  of  pure  water  resting  on  a clean,  flat  surface  of  PTFE.  The  red 
dashed  line  shows  the  profile  calculated  using  the  same  physical  parameters  as  in  the  experiment. 
Photograph  courtesy  of  F.  Barberis  [2]. 

energy  gave  7 iv  = 72  mJ/m2,  with  very  little  change  as  a function  of  time.  Figure  5 shows  a 
photograph  of  a side  view  of  a 2-mL  drop  from  that  study.  The  height  of  the  drop  was  measured  10 
times,  giving  an  average  value  of  4.61  mm  and  an  uncertainty  of  0.2  mm,  expressed  as  one  standard 
deviation  [2].  The  macroscopic  contact  angle  was  measured  with  image  analysis  software  and  was 
reported  to  be  104°  ± 0.5°,  where  the  uncertainty  is  expresssed  in  terms  of  the  precision  of  the 
measuring  device  [2]. 

Superimposed  on  the  photograph  in  Fig.  5 is  a calculated  equilibrium  profile.  For  the  calcula- 
tions, the  conditions  of  the  experiment  were  reproduced  as  closely  as  possible  (7 iv  = 72  mJ/m2, 
p = 1000  kg/m3,  and  g = 9.8  m/s2).  Different  values  of  Qw  were  calculated,  and  the  best  fit  to 
the  experimentally  observed  profile  was  obtained  for  0^  — 110°.  This  is  a modestly  greater  value 
than  that  reported  experimentally  [2].  However,  it  should  be  remembered  that  only  the  apparent, 
macroscopic  contact  angle  was  measured.  The  distortion  of  the  drop  due  to  the  influence  of  gravity 
can  cause  the  apparent  contact  angle  to  deviate  somewhat  from  the  actual  thermodynamic  contact 
angle,  the  latter  of  which  can  reliably  be  measured  only  at  much  higher  magnifications.  Taking 
these  considerations  into  account,  the  agreement  between  theory  and  experiment  is  remarkably 
good. 

Eq.  (60)  indicates  that  the  pressure  within  a liquid  drop  should  be  a sharply  decreasing  function 
of  the  drop  volume.  Fig.  6 shows  this  relation  for  water  (7 iv  = 72.2  mJ/m2)  when  Qw  = 45°,  0^  = 
90°,  and  0^  = 135°.  As  expected,  the  pressure  difference  increases  with  increasing  contact  angle 
because  the  drop  assumes  an  equilibrium  shape  with  increasing  mean  curvature  as  Ow  increases, 
as  shown  in  Fig.  3. 

4 Discussion  and  Conclusions 

Initially,  the  analysis  presented  in  this  report  was  undertaken  with  the  hope  of  comparing  Eq.  (7), 
recently  proposed  by  Searcy  [1],  to  the  more  commonly  accepted  form  in  Eq.  (2).  It  quickly  became 
evident  that  the  only  difference  between  these  two  thermodynamic  potentials,  from  a macroscopic 
point  of  view,  is  the  slightly  different  expression  that  is  derived  for  the  thermodynamic  contact 
angle  given  in  Eq.  (57).  In  fact,  as  already  described,  the  simple  substitution  A/  = (7 is  — 7 sv)  — 7 iv 
makes  that  equation  identical  to  the  equation  of  Young  [9].  Other  than  this,  the  consequences  of  the 
Searcy  formulation  and  the  Gibbs  formulation  are  mathematically  equivalent.  That  is,  there  is  no 
way  to  distinguish  between  the  two  thermodynamic  potentials  using  only  macroscopic  observations 
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Figure  6:  Pressure  difference  across  the  liquid-vapor  surface  as  a function  of  volume  of  sessile  liquid 
drops  when  7 iv  = 72.2  mJ/m2,  p = 1000  kg/1113,  and  g = 9.8  m/s2. 
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of  drop  shapes.  Adamson  [10]  has  pointed  out  that  there  are  no  apparent  experimental  means  to 
measure  7 is  and  7 sv  independently.  In  fact, if  7 iv  is  known  from  a separate  experiment,  measurement 
of  the  contact  angle  in  a wetting  experiment  can  provide  only  the  difference  between  7 /s  and  ysv. 
Searcy’s  A 1 parameter  provides  another  way  to  look  at  the  wetting  phenomenon  that  basically  avoids 
the  issue  of  the  individual  values  of  7 is  and  7 sv . In  principle,  the  thermal  free  energy  contribution 
to  A 1 could  be  independently  measured  by  immersion  calorimetry,  but  there  are  also  likely  to 
be  configurational  entropic  contributions  that  likely  would  make  an  independent  measurement  of 
A 1 difficult.  Nevertheless,  despite  the  conceptual  differences  between  Searcy’s  approach  and  the 
traditional  approach  to  wetting,  this  report  demonstrates  that  they  both  make  precisely  the  same 
predictions  about  the  equilibrium  shapes  of  sessile  drops  on  rigid  solid  substrates. 

The  common  approach,  employed  by  Young  [9],  for  relating  the  contact  angles  is  to  view  jsv, 
7 is,  and  7 iv  as  force  vectors  that  must  balance  at  mechanical  equilibrium.  The  Young-Dupre 
equation  follows  by  summing  the  horizontal  components  of  these  three  forces  to  zero.  In  contrast, 
the  present  analysis  shows  how  Eq.  (1)  arises  as  a necessary  boundary  condition  imposed  by  the 
minimization  of  the  system’s  thermodynamic  potential.  The  same  conclusions  have  been  reached 
using  variational  methods  on  other  systems  [11, 12],  I11  the  usual  force-balance  approach  for  deriving 
the  Young  equation,  the  balance  of  forces  in  the  vertical  direction  is  neglected.  It  seems  to  be 
implicitly  assumed  that  the  solid  substrate  is  capable  of  balancing  the  vertical  component  of  the 
liquid-vapor  surface  tension.  In  the  energy-minimization  approach  used  in  the  present  study,  the 
analogous  assumption  is  that  the  substrate  is  rigid,  i.e.  it  stores  no  elastic  energy  because  it  cannot 
be  deformed.  In  reality,  all  solids  are  deformable,  and  at  equilibrium  the  substrate  surface  will 
adopt  a shape  that  minimizes  a modified  thermodynamic  potential  in  which  the  elastic  energy  of 
the  substrate  is  included  along  with  all  the  other  contributions  that  have  been  considered  in  this 
analysis.  I11  the  limiting  case  for  which  the  substrate  has  zero  shear  modulus  (i.e.  a liquid)  it  has 
been  shown  [12]  that  the  condition  at  the  three-phase  junction  is  that  of  a von  Neumann  triangle 
when  gravity  is  neglected,  regardless  of  whether  the  solution  is  obtained  by  minimization  of  the 
thermodynamic  potential  or  by  a force  balance  approach.  Similar  results  also  have  been  obtained 
for  liquid  drops  on  a thin  elastic  membrane  [13]. 

In  the  present  case,  it  has  been  shown  that  gravity  has  no  influence  on  the  thermodynamic 
contact  angle,  although  clearly  gravity  can  influence  the  macroscopic,  or  apparent,  contact  angle 
observed  at  lower  magnifications.  More  generally,  no  body  forces  can  exert  an  influence  on  the 
thermodynamic  contact  angle.  Previously,  it  also  was  shown  that  the  thermodynamic  contact 
angle  is  not  influenced  by  elastic  mismatch  stresses  between  a wetting  solid  particle  and  its  solid 
substrate  [14, 15].  Therefore,  the  contact  angle  is  a robust  property  of  a given  system,  and  depends 
only  on  the  energy  densities  of  the  participating  interfaces. 
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A Appendix 

In  this  section  we  derive  the  form  for  the  mean  curvature  of  the  surface  shown  in  Fig.  1.  In  spherical 
polar  coordinates,  the  position  vector  x specifying  the  surface  is  parameterized  by  0 and  0.  We 
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will  restrict  attention  to  surfaces  of  revolution  like  that  considered  in  this  study,  for  which 


x (cf),  9) 


/ 7 Z(0)  sin  0 cos  0 \ 
71(6)  sin  9 sin  (f> 

V 71(6)  cos  6 


A.l  First  Fundamental  Form 

The  first  fundamental  form  for  the  arc  length,  s,  is  then  given  by  [6] 

d s2  = g (d</>)2  + 2g(ped(f)d6  + gee  (d6)2 
where  the  components  of  the  tensor  g are 


dx  dx 
d(f>  d(f> 
n2  sin2  6 


dx  dx 


g<t>e  = - 


d(f>  d6 
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dee 


dx  dx 
d6  d6 
7Z2(6)  + 1Z2(6) 


A. 2 Surface  Normal  Vector 

The  unit  vector  h(4>,  6)  normal  to  the  surface  at  any  point  is  given  by 


n = 


dx  v dx_ 
dd  x dd 


dx  dx_ 
dd  x dd 


Evaluating  the  derivatives  in  Ecp  (66)  and  simplifying  gives 


7Z2  + ^ 


1/2 


n cos  0 cos  4>  — TZ  sin  0 cos  d>  \ 
1Z  cos  0 sin  — TZ  sin  0 sin  <p 
V 1Z  sin  0 + 7Z  cos  0 J 


(61) 
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A. 3 Second  Fundamental  Form 

The  second  fundamental  form,  II , relates  to  the  local  deviation  of  the  surface  away  from  its  tangent 
plane  at  the  point  of  tangency  [6]: 


II  = (d0)2  + 2&0fld0d6>  + bee  (d#)2  (68) 

where  the  coefficients  of  the  tensor  b are 


R2  sin2  9 — RR  sin  9 cos  9 
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(n2  + iz2j 


(71) 


A. 4 Mean  Curvature 

The  mean  curvature  of  any  smooth  surface  is  given  in  terms  of  the  tensors  g and  b according  to  [6] 


r 9<p4>Iee  - 2 g^el^e  + 9ee^<t>4> 

K,  = x — 

9cj)(p9ee  ~ 9§q 


(72) 


For  a surface  of  revolution,  this  can  be  simplified  considerably  using  the  results  of  the  previous  two 
sections: 


iz  — iz  cot  o t iz2  + in2  - nit 

77  \ 1/2  77  \ 3/2 

n (n2  + r2)  (n2  + n2j 


(73) 


The  boundary  condition  on  the  surface  at  9 = 0,  given  in  Eq.  (46),  is  that  7t(0)  = 0.  Therefore, 
the  term  containing  cot  6 in  the  previous  equation  will  be  well-behaved  in  the  neighborhood  about 
9 = 0 as  long  as  R approaches  zero  faster  than  tan  9 approaches  zero.  We  can  examine  two  limiting 


18 


cases.  First,  for  small  enough  drops,  experiments  show  that  the  surface  is  shaped  like  a section  of 
a sphere  with  radius  rs.  In  this  limiting  case, 


no) 


rs  (1  — cos  Qw)  sec  9 


i - (i 


cos  sec2 


tan  6 


The  tan  9 in  this  expression  cancels  cot  9.  so  the  first  term  in  Eq.  (73)  remains  finite  as  6 approaches 
zero. 

Another  limiting  case  is  that  of  very  large  drops.  Experimentally  obtained  cross-sectional 
profiles  on  large  drops  indicate  that  they  are  approximately  planar  near  the  center  of  symmetry  (see, 
for  example,  Fig.  5.  If  the  drop  is  planar  in  a neighborhood  near  0 = 0,  then  in  this  neighborhood. 


77(0)  ~ 77  ( 0)  sec  6 


The  the  derivative  with  respect  to  6 is 

77(0)  ~ 77.(0)  sec  0 tan  9 

Again,  the  tan  6 in  this  expression  and  cot  9 terms  cancel,  so  the  first  term  in  Eq.  (73)  remains 
finite  when  the  drop  is  nearly  flat  at  its  apex. 

At  the  point  9 = 0,  the  assumption  of  rotational  symmetry  dictates  that  the  two  principal 
curvatures  must  be  equal.  Therefore,  for  any  drop  size  or  shape,  the  mean  curvature  at  0 = 0 must 
be  given  by 


/C(0)  = 2 


966 

bee 


2 77 2 + 4772  - 27777 


77'2  + 77'2 


3/2 


(74) 


B Appendix 

Because  cot  9 in  Eq.  (59)  is  singular  at  9 = 0.  a change  both  in  parameterization  of  the  liquid-vapor 
surface  and  a transformation  of  the  coordinate  system  is  used,  as  shown  in  Fig.  7. 

The  relations  between  quantities  in  the  old  and  new  coordinate  system  are: 


sin  <p  — 

- nt ) 

(75) 

no)  - 

- -m 

(76) 

©- 

- ntf) 

(77) 

no)  - 

- -ntf) 

(78) 

where  t = s/l  is  the  dimensionless  arc  length,  tf  is  the  total  arc  length  of  the  drop  profile,  and 
where  the  new  functions  77(s)  = r(s)/l  and  77 (s)  — z(s)/l  are  understood  to  be  dimensionless. 

It  is  readily  shown  [16]  that,  in  the  new  coordinate  system,  the  principle  curvatures  are 
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+z 


Figure  7:  Cylindrical  polar  coordinate  system  used  in  the  numerical  solution  of  the  Euler-Lagrange 
equation.  The  liquid-vapor  surface  is  parameterized  by  its  arc  length  s.  ip(s)  is  the  tangent  angle 
to  the  drop  (ip  < 0). 


1 dip 

R\  ds 

1 sin  ip 

~R2  = _ d R 


(79) 

(80) 
(81) 


Because  the  mean  curvature  is  invariant  under  coordinate  transformations,  we  may  substitute  this 
result  into  Eq.  (59).  Thus  the  Euler-Lagrange  equation  (59)  is  transformed  to 


V = (3Z(t) 


dip 

~dt 


sin  ip 

n 


(82) 


In  this  equation,  we  use  rj  = 77  + (3  Z ( tf ).  Furthermore,  from  Fig.  7,  the  tangent  angle  ip  is  related 
to  Z and  7Z  according  to 


dIZ 

d T 
dZ 

~dt 


cos  ip 
sin  ip 


(83) 

(84) 

(85) 


Eqs.  (82)-(84)  thus  are  a system  of  three  coupled  first-order  ordinary  differential  equations  that 
must  be  solved  simultaneously  for  Z(t),  TZ(t),  and  ip(t).  Furthermore,  using  Eqs.  (75)-(78),  the 
boundary  conditions  given  in  Eq.  (46)  and  (56)  are  transformed  to: 
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0 (tf)  = -Qw 

(86) 

72.(0)  = 0 

(87) 

Z{  0)  = 0 

(88) 

(89) 

In  Eq.  (82),  the  term  sin  0/72.  has  a singularity  at  t = 0,  which  can  be  removed  by  expanding  72. , 
Z,  and  0 about  that  point: 


m)  = rit  + p/ + ’f/ + + 0(t9)  (90) 

z(f)  = li2  + ii4  + ^6  + 0(t8)  (91) 

m = i’lt  + ~ti  + ^t5  + ^t7 + 0(t9)  (92) 

(93) 

which  assumes  that  7 Z(t)  and  00)  both  are  odd  functions  of  t,  and  that  Z(t)  is  even.  The  coefficients 
in  these  expansions  are  determined  by  substituting  the  expansions  into  Eqs.  (82),  (83),  and  (84). 
From  Eq.  (83),  we  find 


n = 1;  r3  = 

r5  = -4030i  + ip\ ; r7  = -60 50i  - 1O03  + 20 0303  - 0^ 


From  Ecp  (84), 


32  = '00 


^4  = 03  - 0i ; 


^6  = 05  - 100301  + 0] 


From  Eq.  (82), 


01  = 


zR. 

2 ’ 


03  = 


-3/7  f) 


05  = 


(3f  5/32fj 


-11(3  f 33(32f  35/73t) 

07  = + 


16  16  ’ r‘  256  64 

Substituting  these  coefficients,  the  singular  term  sin 0/72.  becomes 


128 


sin0  —fj  (3fjt 


n 


— 1 7~  + 


2 


16 


(3f-(3zil 

384 


0 + 


2-3 


43o 


61440  46080  18432 


PVS 5 +31/3V7 


t6  + 0(18) 


(94) 


This  expression  is  substituted  for  sin  0/72.  over  a small  interval  0 < t < to  in  defining  the  ODE  for 
numerical  solution.  If  to  = 0.01,  the  expansion  should  give  14  or  15  digits  of  accuracy. 
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